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Sound  Propagation  in  the  Sea:  — 

THE  FINITE  CHANGE  IN  WAVEFRONT. DIVERGENCE  IN  PASSING  THROUGH 
"discontinuities  IN  sound-speed  GRADIENT  AS  EFFECTED  AT 
REFLECTIONS . 


Sound  Propagation  in  the  Sea: 

THE  TINITE  CHANGE  IN  WAVEFRONT  DIVERGENCE  IN  PASSING  THROUGH 
DISCONTINUITIES  IN  SOUND-SPEED  GRADIENT  AS  EFFECTED  AT 
REFLECTIONS 


1.  Definitions 


We  focus  our  'httention  on  the  problem  of  sound  intensity  mapping  in 
a vertical  plane. ^ Tor  a point  source,  each  selected  ray  is  associated  with 
an  angle,  V , , at  which  the  ray  emanates  from  the  source;  by  convention 
this  angle  is  measured  clockwise  from  the  horizontal.  The  labelling  of  rays 
by  their  source-emission  angles,  , defines  y wherever  the  medium  is 

u b 

traversed  by  the  rays.  It  should  be  noted  that  is  not  necessarily  single 
valued  throughout  the  medium;  more  than  one  ray  may  pass  through  a point. 

Such  multiple  values  are  associated  with  foldings  of  the  wavefront.  We  may  If 

refer  to  a spread  of  rays  extending  from  one  fold  to  another  as  a "family"  of 

rays. 

For  a family  of  rays,  in  the  vertical  plane,  the  specific  wavefront 
length,  denoted  by  L , is  defined  at  each  point  by 


where  « is  a unit  vector,  clockwise  normal  to  the  ray  direction  at  that  point. 

» i 

Negative  L thus  Implies  V is  decreasing  in  the  direction  m;  sign  changes 
in  L are  important  indications  of  foldings  of  the  wavefront  as  at  caustics  or 
by  reflections. 

• » 

The  differential  equation  for  integrating  L along  a ray  is  developed 
in  Reference  [l]  as  , 
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where  s Is  ray  path,  C Is  sound  spued  and  n Is  ray  normal.  The  Initial 
conditions  at  the  source  point  are 


Le  - 0 , (diyds)c  « 1 
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The  patametei  1.  may  be  integrated  .iloiui  each  tay  beginning  at  the 
source.  The  advantages  ot  such  diiect  integration  include  (l)  the  pie- 
cluslon  ot  the  problem  ot  measuring  specific  spreading  by  geometnc 
measurement  ovei  t i i\i t distances  between  dlver.se  lays,  (2)  the  iudicat ton , 
at  lorded  by  sign  changes  in  l,,  ot  foldings  ot  the  wavetiont  caused  by 
reflections  oi  at  caustics,  and  (a)  the  indication,  attoided  by  the  spatial 
continuity  in  l.,  of  the  adequacy  ot  the  intensity  resolution  attoided  by  the 
selected  family  ot  lays. 

Turning  to  t'tq.  1 we  note  that  the  ret  lectori  lays,  attei  suitace 
reflection  at  A,  have  increasing  to  the  left  of  the  ray  direction, 
implyinq  negative  I.  by  convention.  To  accomodate  the  convention  the 
sign  of  1.  must  be  changed  at  the  point  of  icfleetion , by  special  llioyision 
in  the  integration  of  Eq,  (2)  along  a ray.  At  B,  1.  passes  through  .-i'io, 
again  changing  sign;  this  sign  change  is  inhei out  in  the  integiation  ot 
Eg.  (2)  along  a ray.  Such  foldings  of  the  wavetiont  in  the  medium  are 
called  caustics;  the  caustic  Is  the  locus  of  L 0 from  ray  to  lay. 

Not  only  is  special  provision  inquired  .it  reflection  to  change  the 
sign  of  L but,  generally,  also  the  spieading  tale  requites  special  tieat- 
nient.  At  reflection  the  results  are  such  as  they  would  he  it  the  medium 
weie  to  continue  in  mlrroi  image  beyond  the  reflecting  surf  act',  as  shown 
by  the  dashed  sound-speed  profile  and  the  dashed  continuation  of  the  lays. 


For  the  sound  speed  having  a gradient  normal  to  the  reflecting  surface,  the 
ray  effectively  passes  through  a kink  (gradient  discontinuity)  in  the  profile 
at  the  point  of  reflection. 

The  integration  of  Eq.  (2)  along  a ray  can  only  accomodate  the  effect 
of  continuous  gradients  of  the  sound  speod,  C.  Thus  special  treatment  is 
necessary  in  passing  through  gradient  discontinuities;  the  spreading  rate, 
dL/c)s  , abruptly  changes.  We  shall  develop  the  treatment  for  gradient 
discontinuities  — interfaces  within  the  medium  — and  adapt  the  treatment 
to  reflection. 


Fig.  1 The  traces  of  two  rays,  “X;  ■ 0°  and  ?g  ■ 5°,  leaving  a point 
source,  in  a sound  speed  distribution  which  increases  linearly 
with  depth. 


2.  Treatment  at  Interfaces 


Depending  on  the  scale  of  such  features  of  — and  on  the  resolution 
of  the  numerical  representation  of -- the  sound  speed  distribution,  it  may 
be  warranted  and  even  desirable  to  treat  strong  changes  in  gradient  by  an 
interface  of  discontinuity. 


fig.  2 An  interface  (dashed)  of  discontinuity  in  sound-speed  gradient 
depicted  by  speed  isolines. 


In  Fig.  2 we  depict  an  interface  of  discontinuity  in  sound  speed 
gradient.  In  the  seas,  the  slope  of  such  interfaces  is  generally  negligible. 
In  our  development  which  follows  we  shall  assume  the  ray  to  meet  the 
interface  at  depth  Z*  and  we  neglect  the  slope  of  the  interface. 

In  Fig.  3A  we  depict  a ray,  crossing  the  interface  at  depth  Z*  , 
with  an  arbitrary  angle,  , to  the  horizontal.  In  order  to  derive  the  effect 
we  expand  the  interface  into  a transition  zone  as  shown  in  Fig.  3B.  The 
effect  may  be  determined  by  integrating  Fq.  (2)  through  the  transition  zone, 
from  Z*  to  ZA  + h,  and  then  letting  h tend  to  zero. 


I NTERFACE 
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Fig.  3 Discontinuity  interface  at  Z*  (Fig.  3A) , expanded  into  infinitesimal 
transition  zone,  with  continuous  gradient,  of  depth  increment  h 
(Fig.  3B) . 


The  sound-speed  profile,  in  the  transition  zone,  is  fitted  by 


c = c*  + 
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az 
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(4) 


which  accomodates  the  value  C*  at  Z^  . Differentiation  yields  the  gradient: 
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which  is  continuous,  going  from  the  value  just  before  the  interface, 
(ac/az)_  , to  that  just  after  the  interface,  (aC/az)+  . 

We  express  the  integration  of  Eg.  (2)  through  the  transition  zone  by 


s*  + h/sin  y* 
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Since  we  shall  let  h go  to  zero  we  anticipate  the  simplifications  and 
reduce  Eq.  (6)  to 


s*  + h/ sin  y* 
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ds 
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For  the  transition  zone  we  have 
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which  is  independent  of  s.  Thus  Eq.  (7)  yields 
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which  is  unchanged  as  h -*  0. 


3.  Formula  for  Reflections 

For  reflections  we  have  the  added  complication  that  the  sign  of  L 
is  flipped: 


L 


R 


(10) 


The  subscript  I is  used  to  denote  an  incident  value  and  the  subscript  R a 
reflected  value.  Taking  this  into  account,  the  spreading  rate  changes 
according  to 
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For  numerical  treatment  we  substitute 


(12) 
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LR  ~ LR  - 1 
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and  Eq.  (10)  into  Eq.  (11)  to  obtain 
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The  new  values,  LR  and  LR_1  , are  used  in  place  of  Lj  and  in 

proceeding  with  the  integration  of  Eq.  (2)  along  the  ray  following  reflection. 
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